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Abstract Allozymes of eleven species of Megophryinae in China were examined 
electrophoretically to investigate genetic diversity and phylogenetic relationships. Fourteen 
enzymes, presumptively coded by 24 loci were detected to be variable. Gene frequencies of 
each population at each locus were presented. The commonly used measure of genetic 
diversity, the average heterozygosity (H) were calculated based on gene frequencies. The re- 
sults indicated that Megophryinae had a high level of genetic diversity in amphibians, an aver- 
age H of 0.18, ranging from 0.058 to 0.28. Nei’s (1978) genetic distances(Nei? s D) were calcu- 
lated for all possible population pairs. A dendrogram of 13 populations representing 11 
species, 3 genera of Megophryinae were derived and presented by using UPGMA, based on 
Nei s D. The assignment of Ophryophryne as a distinct genus were supported by an average 
Nei! s D of 1.4067 which separated O. microstoma from all other populations.Subdivision of 
Brachytarsophrys from Megophrys was not supported by this study. Within Megophrys, 
three groups were recognized: (D M. lateralis, M. giganticus and M. longipes;@ M. 
palpebralespineosa, M. boettgeri and M. parv; @ M. minor and M.  kuatunensis. 
Three populations of M. omeimontis were closely related and share a clade independent from 
all other Megophrys, and B. feae as well. 


Key words Amphibia, Pelobatidae, Megophryinae, Genetic diversity, Phylogenetic relation- - 
ships, Allozyme electrophoresis, China . 


1 Introduction 


Genetic variation is one of the fundamental parameters of evolutionary process. 
This is because the potential of a population is a function of the amount of genetic varia- 
tion present in the population at a given time. The molecular approach has made possi- 
ble an extensive comparison of organisms at levels that closely reflect gene products (Des- 
sauer, 1974). Biochemical methods of both protein electrophoresis and DNA 
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sequencing data has been used with increasing frequency to taxonomic and evolutionary 
problems in amphibians (see review of Nevo and Beiles, 1991). Starch gel 
electrophoresis (SGE) is a fairly common and simple technique that makes possible the 
study of protein variation with only a moderate investment. In this study, we present a 
allozyme electrophoresis approach to 11 species of Megophryinae in China. 

The oriental megophryines are much diverse in external appearance, they are terres- 
trial or semifossorial and some have fleshly dorsolateral folds and / or supraeiliary pro- 
cesses. The Chinese megophryines were subdivided into two subfamilies by Dubois 
(1983, see Frost, 1985), Tian et al (1983). The subfamily Megophryinae governed 
four genera, i. e. Megophrys, Brachytarsophrys, Atympanophrys and Ophryophryne. 
Among them, the genus Megophrys is highly diversified, with over twenty species dis- 
tributed in Southeastern Asia from India to China, south to Philippines, Sumatara, 
Java, Borneo and Natuna. Over half of its members occurred in China, some of them 
are endemic to China. The taxonomic and systematic problems of this group are far 
from being worked out. For example, the genus Ophryophryne  (Boulenger, 1883) 
were considered synonym of Megophrys by some authors (Frost, 1985; Duellman et 
al, 1986), but considered independent genus by others(Tian et al., 1985; Zhao et 
al, 1993). Further more, Liu et a/.(1961) pointed out that morphological differences 
were existed among geographic populations of M. omeimontis, as well as in same case 
of Brachytarsophrys carinensis and M. minor. Hu (personal discussion, 1992), after 
comparing external and internal morphology of B. feae and M. montana (type spe- 
cies of Megophrys), suggested that they resembled each other. 

The proposal of this work is not to solve all the above problems which may need 
considerable many works. This study is carried out to begin this process of under- 
standing. Two questions are addressed: (D What is the phylogenetic relationships within 
this group? 29 How much does the genetic diversity exist in the members of this group? 


2 Materials and Methods 


Between May, 1991 and August, 1993, a total of 188 specimens representing 13 
populations of Megophryinae was collected during field trips; information on localities 
and sample sizes is given in table 1. Tissues of heart, liver and skeletal muscle were cut 
out immediately after sacrificing the animals and put in the liquid nitrogen(-196'C), and 
then brought to laboratory and maintained at —70C in refrigerator until use. Voucher 
specimens are deposited in Collection of Kunming Institute of Zoology, the Chinese 
Academy of Sciences. 

Extracts were prepared by homogenizing each sample after method presented by 
Pasteur et al.(1988), and maintained at -10C until use. Fifteen specific enzymes were 
assayed by horizontal starch gel electrophoresis (SGE), and one, i. e. «-amylase, 
was assayed by polyacrylamide gel electrophoresis PAGE). Enzymes, electrophoretic 
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conditions and staining references are listed in Table 2. 


Table 1 Localities and sample sizes of thirteen populations of Megophryinae 





` Population and species 


1 M. 
2 M. 
3 M. 
4 M. 
5 M. 
6 M. 


lateralis 
omeimontis 
omeimontis 
omeimont is 
minor 


kuatunensis 


7 M. 
.8 M. 
9 M. parva 
10 M. longipes 
11 M. gigantica 
12. O. microstoma 
13 B. feae 


palpebrales pineosa 
boettgeri 


Locality 


Tengchong, * Yunnan 
Mt. Emei, Sichuan 
Jingdong, Yunnan 
Pinbian, Yunnan 
Jingdong, Yunnan 
Dehua, Fujian 


| Chun-an,Fujian 


Pinbian, Yunnan 
Chun-an, Yunnan 
Mengla, Yunnan 
Hekou, Yunnan 
Jingdong, Yunnan 
Hekou, Yunnan 
Jingdong, Yunnan 


Elevation MS 
i Sample size 
(approx. in m) 
2340 20 
750 | 20 
2000 20 
700—1100 3 
2000 . — 20 
1300 2 
1180 3 
700-1100 - | 
1180 15 
600- 700 12 
120- 200 10 
2000 20 
120- 200 4 
2000 22 


Table 2 Enzymes, electrophoresis conditions and staining references 


Enzyme 
Acid phosphatase 
Alcohol dehydrogenase 
Estrase 
Glucose dehydrogenase 
Glutamate dehydrogenase 
L-Lactate dehydrogenase 
Malate dehydrogenase 
Alkaline phosphatase 
a-Amylase 
Octanol dehydrogenase 


Peroxidase 


Catalase 
Sorbitol dehydrogenase 


Superoxide dismutase ' 


E. C. 
Number 
3.1.3.2 


1.1.1.1 
3.1.1.1 
3.1.1.2 
1.1.1.47 
1.4.1.2 
1.1.1.27 
11.137. 
3.1.3.1 


3.2.1.1 
1.1.1.73 


1.11.1.6 
1.1.1.14 


1.15.1.1 


Presumptive 


locus 
ACP-1 
ACP-2 
ADH-1 


ADH-2 


EST-1 
EST-2 
GLC-1 
GLC-2 
GLD-1 
GLO-2 
LDH-1 
LDH-2 
MDH-1 
MDH-2 
AKP-2 
AMY-1 
ODH-1 
ODH-2 
PER-1 
PER-2 
CAT-1 
SDH-1 
SDH-2 
SOD 


Electrophoresis type Staining 
and condition reference 
SGE,Tris-Citrate p''* 
pH7.0 

SGE,Poulik " | P 
SGE, Poulik © P 
SGE,Tris-Citrate P 
pH7.0 

SGE, Tris—Citrate P 
pH7.0 

SGE,Poulik and P 
Tris-Cirtrate,pH 7.0 

SGE, Tris-Citrate P 
pH7.0 

SGE, Tris-Citrate, pH7.0 P 
PAGE, Tris-Gly' " P 
SGE, Poulik P 
SGE, Tris-Citrate, S&P'''" 
pH7.0 

SGE, Poulik S&P 
SGE, Poulik P 
SGE, Poulik P 





* Poulik- Tris-Citrate-Borate, pH8.7; * * discontinues, pH6.7 and 8.9; 


+ + * Pasteur et al., 1988; 


* * * * Shaw and Prasad, 1970. 
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Because malic enzyme(NADP-dependent malate dehydrogenase) and aldenyde 
oxidase were not considerably interpretable, these two enzymes were discounted from 


analysis. Thus, a maximum of 14 enzymes coded by 24 loci were used. : 





Figure 1 Some of the photographs of gels obtained by isozyme electrophoresis 


Il. AMY(a-Amylas) 2. LDH(L-Lactate dehydrogenase); 3. PER (Peroxidae); 
4. EST(Esterase); 5. AKP(Alkaline phosphatase); 6. GLD(Glutamate dehydrogenase): 
7. GLC(Glucose dehydrogenase); 8. ADH(Alcohole dehydrogenase) 
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Table 3 Allelic frequencies for variabe loci 
Populations 
.. Locus 1 2 3 4 5 6 7 8 9 0 11 12 B 
ACP-1 A 0.50 0.44 0.33 1.00 | 
B 0.85 0.50 056 067 100 — 0.95 1.00 
C 015 1.00 0.05 0.10 0.83 1.00 
D © «0990 017 
E I | 1.00 
ACP-2 A 1.0 1.00 923 «7 | (c 
| 072 T" ... 9,65 0.90 1.00 1.00 
C 0.11. L00 1.00 1.00 |. 100 012 0.10 
D 017 , E T" 1.00 
 ADH-1A 100 042 0.60 0.67 -© 0.47 0.15 0.37 0.56 
| B 0.58 1.00 1.00 0.40 0.33 1.00 0.3 085 | | 0.44 
C | 1.00 0.63 
D . P E 1.00 
ADH-2A 1.00 0.38 1.00 1.00 0.44 0.05 0.93 0.30 1.00 
B 0.62 100 | 100 100 0.566 0.95 * 0.07 0.70 1.00 
AKP-1 A 0.25 030 030 © © 633 023. 0.05 
B 0.75 0.70 0.70 0.38 . 0.56 0.67 0.64 0.50 0.10 1.00 
C 100 0.62 1.00 0.44 0.13 0.50. 0.85 1.00 
AMY-1 A 1.00 
B 1.00 1.00 1.00 1.00 1.00 0.10 0.08 
C 087 0.90 0.92 0.95 017 037 
D 013 E | 0.05 0.83 0.63 
E : 1.00 
EST-1A 046 ^ ' 6.65 0.56 0.80 0.45 0.50 0.10 1.00 021 
B 0.4 1.00 1.00 1.00 1.00 035 044 0.20 0.55 0.50 0.90 | 0.79 
EST-2A | ! 1.00 
B 100 0.67 0.95 0.46 0.70 0.60 0.50 0.47 1.00 
C 1.00 0.33 0.05 0.54 0.30 0.40 0.50 0.53 
D 1.00 1.00 
GLC-1A 0.15 0.83 |. 025 020 0.79 0.54 0.83 
B 0.58 0.17 100 1.00 0.75 0.70 .021 085 0.46 0.17 0.71 0.50 
C 0.27 010 — 0.15 |. 0.29 1.00 0.50 
GLC-2A 013 046 100. 008 100 100 100 0.80 1.00 
| B 0.87 0.54 100 092 ^. | 1.00 1.00 020. 0.50 
| € Ar 8 P | 0.50 
GLD-1.A 036 035 0.5. 0.83 017 038 0.56 0.60 0.36 0.60 0.50 0.19 
B 0.64 0.54 0.10 0.17 0:78 0.50 0.44 0.7 0.20 0.64 0.40 0.50 0.81 
C 0.11 0.85 0.09 0.12 0.23 0.80 
GLD-2A 100 100 1.00 1.00 1.00 0.50 
B | ^» 1,00 1.00 1.00 1.00 1.00 1.00 0.67 0.50 
C 0.33 
LDH-1A 0.05 | 1.00 
B. 095 1.00 0.08 1.00 1.00 1.00 0.19 1.00 
C 1.00 1.00 1.00 0.92 0.81 | 
D 


1.00 
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(continued) 
Populations 
~ Locus 1 2 3 4 5 6 7 8 9 0 mH 2 B 
LDH-2 A 0.10 s 1.00 | 
B 1.00 069 1.00 1.00 090 022 1.00 1.00 1.00 1.00 1.00 ' 1.00 
C 4 031 . 0.78 
MDH-1A 040 1.00 0.10 | | 0.30 0.44 0.18 0.67 
B 0.60 031 1.00 0.60 1.00 1.00 0.82 0.70 0.56 0.61 0.33 0.29 
C 0.69 0.30 . 0418 0.21 0.71 
MDH-2 A 0.0 025 | 0.60 0.70 1.00 9.11 0.07 0.67 
B 100 0.27 0.12 0.50 040 0.30 1.00 1.00 0.89 0.72 0.33 1.00 
C 0.11 0.63 0.50 0.21 
PER-1A 1.00 1.00 1.00 1.00 1.00 0.10 1.00 0.05 1.00 0.33 
B | 0.90 1.00 0.95 0.67 
C 1.00 1.00 
PER-2 A 1.00 
B 073 | 0.25 1.00 0.78 1.00 0.75 
C 027 100 1.00 1.00 | 0.85 0.43 0.11 0.25 
D 1.00 1.00 0.15 0.32 0.11 
SDH-1A 0.45 1.00 1.00 0.25 0.37 0.53 0.90 0.79 0.68 
B 055 1.00 1.00 0.75 0.63 0.47 0.10 021 0.32 1.00 
CC | | | | 1.00 
SDH-2A 0.55 0.54 1.00 1.00 0.39 0.31 023 0.05 0.50 = 0.10 
B 045 0.46 1.00 0.61 0.69 0.77 0.95 1.00 0.50 1.00 0.90 
ODH-1A 013 0.63 0.38 
B 087 0.50 0.10 0.83 1.00 0.62 1.00 1.00 1.00 1.00 
C 0.37 0.50 1.00 0.90 0.17 1.00 | 
ODH-2 A 0.33 0.50 0.50 0.17 0.39 0.15 1.00 1.00 
B 100 0.67 0.50 050 083 0.61 1.00 0.85 1.00 1.00 
C 1.00 
-SOD A 0.80 0.95 
B 0.05 0.10 1.00 0.20 0.05 1.00 095 0.75 0.75 1.00 
C 1.00 0.95 0.90 0.05 0.25 0.25 0.95 ` 0.75 
D 0.05 0.25 
CATA 0.85 0.40 0.50 
B 0.15 0.50 0.50 0.50 0.44 0.60 
C 1.00 1.00 1.00 0.60 0.50 0.50 0.50 0.56 1.00 0.38 
D 0.40 1.00 0.12 





Allele frequencies were calculated, and a chi-square test were used to test for devia- 
tion from Hardy—Weinburg equilibrium. The Nei’s genetic distance (D) was calculated: 
for each possible population pair. Because sample sizes were relatively small (most were 
between 10 and 20), the Nei’ s D were calculated by using the modified method for 
small size presented by Nei(1978). 


3 Results 


In most cases, the pattern of bands observed electrophoretically at any given locus 
was similar to that seen in other vertebrates(for some of the gels, see Fig. 1). The gels 
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were interpreted according to Pasteur et al. (1988) and Harris et al. (1976). Gene fre- 
quencies of each population were calculated and listed in the Table 3. Observed 
heterozygosity of each population at each locus was tested, and in no case was the val- 
ue of heterozygosity observed per-locus significantly different from that expected under 
Hardy-Weinburg law. The average heterozygosity (H) of each population were calcu- 
lated, and the result indicated that the level of H of Megophryinae was very high, with 
a mean value of 0.18, ranging from 0.058 to 0.28(see Tab. 4). 
Table 4 Average heterozygosity (H) of 13 populations representing 11 species of Megophryinae 


population 1 2 3 4 5 6 7 8 9 10 11 12 13 
H 0.213 0.253 0.120 0.089 0.117 0.194 0.172 0.239 0.195 0.214 0.285 0.058 0.209 


Nei’ s (1978) genetic distance (D) was calculated for each possible population pair, 
and given in Table 5. A dendrogram was derived from this D matrix by using 
UPGMA(Fig. 2). We can see that Ophryophryne microstoma share one clade, sepa- 
rate from other populations with an average genetic distance of 1.4067; the remaining 
populations share the other clade, within this clade, the average genetic distance be- 
tween any population pair is below 1.0. "Therefore, this result supports the subdivision 
of Ophryophryne from Megophrys as a distinct genus. But it is not the same case for 
Brchytarsophrys. | 


Table 5 Nei's genetic distance of 13 populations representing 11 species of Megophryinae 
1 2 3 4 5 6 7 8 9 10 11 12 13 


l .0 0.8300 0.8163 0.8460 0.7413 0.9272 0.4998 0.5480 0.6534 0.5253 0.3145 0.1330 0.4892 
2 0 0.3046 0.3573 0.6689 0.8842 0.6662 0.7944 1.1101 0.9857 0.9091 1.3037 0.8038 
3 0 03125 0.8127 0.8226 0.6014 0.6811 0.8003 1.0861 0.9297 1.3602 1.0708 
4 0 0.6913 0.9136 0.7844 0.8860: 0.8000 1.0190 1.1123 1.2825 1.2595 
5 O 0.3384 0.7080 0.4677 0.5299 0.7711 0.6369 1.7702 0.7502 
6 0 0.7844 0.6034 0.9983 0.9306 0.7494 1.5729 0.7573 
7 0 0.3320 0.4386 0.5869 0.6482 0.9906 0.7424 
8 0 .0.4565 0.5125 0.3772 1.1571 0.8097 
9 0 0.5149 0.6892 1.4729 0.9128 
10 0 0.4379 1.3834 0.6403 
11 | 0 1.4447 0.5269 
12 | 0 1.3747 
13 0 


4 | Discussion 


4.1 Heterozygosity 

The average heterozygosity(H) is commonly used as a measure for genetic diversity. 
The levels of genetic diversity measured by H in amphibians are varied non-randomly on 
a massive scale among taxa (families, genera, species), as well as among biotic 
factors( ecological, demographic, and life history). H is structured both taxonomically 
and biotically. H is also higher in amphibians as compared with other vertebra(Nevo et 
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al, 1991). The majority of amphibians hither-to studied had a level of H below 0.13, 
but actually a few taxa had a level of exceeding 0.20 (Green et al., 1988). In the present 
study, H of Megophryinae is very high. This is interesting but difficult to understand, 
. because references of the ecology of Megophryinae is rare. However, Liu et a/(1961) 
pointed out that external morphological divergences were exited in some species of 
Megophrys. And a recent study revealed that skeletal divergences were also existed in in- 
ter- and intra—populations of M. lateralis. So, genetic diversity may be congruent with- 
morphological divergence in Megophrys. 


1 M lateralis 

ll M gigantica | 

10 KM longipes 

1 TT 


8 M boettgeri 


9 M parva 





3 B. /eae 


2.1093 


5 M minor 
6 M kuatunensis 


2 M omeimontis 


3 M ometmontis 


4 M omeimontis 


, | 0.3349 
[RE 12 0 microstoma 


1. 4067 


一 Nei's D 
23.0 1.8 1.6 1. 4 1. 2 1. 0 08 0.6 O04 0. 2 0 


Figure 2 Dendrogram of 13 populations of 11 species of Megophryinae in China, 
based on Nei’s genetic distance, derived by UPGMA 


4.2 Biochemical systematics | 

The family Pelobatidae is formerly divided into two subfamilies, Pelobatinae and 
Megophryinae. Recently, the latter was further arranged into two subfamilies, i.e. 
Megophryinae and Leptobrachiinae(Dubois, 1983, see review of Duellmam, 1986), 
the present study concentrated on Chinese members of the former. The systematics of 
this group are confusing, and the phylogenetic relationships within it is far from being 
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worked out. O. phryophryne was regarded as a distinct genus by Liu et a(1961), Tian 
et al(1985), Fei et a(1990) and Zhao et aK1993), but regarded as subgenus following 
Dubois. (1980). Duellman and Trueb(1986) and Frost(1985) did not subdivided 
Ophryophryne from Megophrys. . Brachytarsophrys carinensis was also divided from 
Megophrys based on external and internal morphological differences from other mem- 
bers of Megophrys(Tian.et al., 1986). But recently, Hu Qixiong( 1992, personal discus- 
sion) after comparing morphology of B. feae with M. montana, suggested that these 
two species resembled each other, and suspected the validity of Brachytarsophrys. 

Thorpe(1983), using Nei’s genetic identity, I, found that members of the same ver- 
tebrate species almost always have identities 0.900 (equivalent to Nei’? s D « 0.105), 
while members of different species within genera characteristically have identities > 0.350 
(equivalent to Nei's D«1.050). (Green et al, 1988). Therefore, the average Nei's D 
of 1.4067 which separate O. microstoma from the remaining species supported the 
subdivision of Ophryophryne from Megophrys asa distinct genus. This result was also 
convinced by tadpole oral characteristics presented by Huang et aK(1991). As for the 
Brachytarsophrys, it is not the same case, an average Nei’s D of below 0.8896 which 
separate feae from other Megophrys did not support the subdivision of feae from 
Megophrys. | 
.. In Megophrys, three groups with intra-group close relationships are recognized: 

(D M. lateralis, M. giganticus and M. longipes, where the former two are closely re- 
lated, with an average Nei! s D of 0.3145, while the latter was separated from the for- 
mer two by an average Nei! s D of 0.4816 © M. palpebralespineosa, M. boettgeri 
and M. parva, where the former two are each other’ s close relatives, with a Nei’ s D 
of 0.3320, an average Nei! s D of 0.4467 separate the latter one from the former two; 

@ M. minor and M. kuatunensis are closely related, with a Nei? s D of 0.3384. 
The above groupings are, at least to some extent, congruent with morphological data 
(Yang et al., in preparation). 

It is special for the case of M. omeimontis. Three geographical populations of this 
species examined here share a clade, this clade separated it from all species of 
Megophrys, and B. feae as well, with a considerable high average Nei! s D, 0.8896. 
Among the three populations of this species, samples from Jingdong, Yunnan and Mt. 
Emei, Sichuan (type locality) are more closely related each other than to samples from 
Pinbian, Yunnan. Because M.  omeimontis is widely distributed in the southwestern 
China, divergences between populations are well known( Liu et al., 1961). Samples | 
from Jingdong, Yunnan are regarded as subspecies by Fei et a/(1990). A further study 
concentrating on various geographic populations of this species is necessary for the fu- 
ture to investigate population divergences and phylogeny. 
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摘要 本文 用 同 工 酶 电 六 的 方法 分 析 了 中 国产 角 蛤 亚 科 3 属 11 种 ， 共 13 个 地 理 种 群 样 
本， 以 探讨 它们 的 遗传 多 样 性 及 系统 关系 。 用 于 分 析 的 酶 共 14 种 ， 由 24 个 基因 位 点 编码 。 
经 过 酶 带 的 分 析 和 计算 ， 得 出 了 每 个 种 群 样本 在 每 个 基因 位 点 上 的 等 位 基因 频率 ( 见 表 3)。 根 
据 基因 频率 的 分 布 ， 计 算得 出 了 各 种 群 样本 的 基因 多 样 性 指数 平均 杂 合 率 (H)( 见 表 4)， 结 
RREH, AAEH H 值 平均 值 为 0.18， 在 两 栖 动 物 中 是 属于 很 高 的 遗传 多 样 性 水 平 ， 根 据 
基因 频率 的 分 布 ， 计 算 了 所 有 种 群 样本 间 的 Ne s 遗传 距 ( 见 表 5)， 并 通过 类 平均 法 
(UPGMA)& @ T 4t Bl (dendrogram). 47 9r 48 JR. 3c Es HW FA M BRL CO phryo phryne)AA fh 8 JR 
(Megoppry9) 中 分 出 为 独立 属 的 观点 ， 但 是 对 短 腿 蟾 属 (Brachytarsophrys) 是 否 为 独立 属 提出 了 
疑问 。 在 角 蛤 属 (Megophrys) 中 ， 可 以 明显 地 看 出 3 个 组 ， 它 们 中 的 种 间 亲 缘 关 系 比 较 近 ， 这 
3 个 组 分 别 是 :中 ELA AR WE CM. lateralis), K tt faM. giganticus)fll K ik fa 9 
(M. longipes) © *H Hz f8&(M. — palpebralespinneosa). . JH FAY (M. boettger) fn [U] Th fA M 
(M. parva @ ;fB8S8(M. minor) FNS fS EOM. kuatunensis), Ria AM. omeimontis)3 
SRR AX, HAAA- ERRER, RSP MI A 25 5838 — 9C. 
我 们 认为 ， 该 种 不 同 地 理 种 群 间 的 比较 还 有 必要 进一步 深入 研究 。 
Xia) PERH, AMV, MeSH, ASL, ICH 





